On the Theory of Metal-Insulator Transitions in Gated Semiconductors 
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It is shown that recent experiments indicating a metal-insulator transition in 2D electron systems 
can be interpreted in terms of a simple model, in which the resistivity is controlled by scattering at 
charged hole traps located in the oxide layer. The gate voltage changes the number of charged traps 
which results in a sharp change in the resistivity. The observed exponential temperature dependence 
of the resistivity in the metallic phase of the transition follows from the temperature dependence 
of the trap occupation number. The model naturally describes the experimentally observed scaling 
properties of the transition and effects of magnetic and electric fields. 



Recently, a metal-insulator transition has been ob- 
served in low-density two-dimensional electronic systems 
- first in Si MOS structures (l]-||], and later in other het- 
erostructures It has been found that when the 

density of 2D electrons n s is below some critical value 
n c s , cooling causes an increase of the resistivity p, while 
at n s > n c s the resistivity decreases with temperature T, 
i.e., the system exhibits an unexpected metallic behavior. 
The insulating phase has been found to be rather usual 
and easy to describe in terms of variable-range hopping 
JlO| . On the contrary, the metallic phase is anomalous in 
at least three respects: i) the p(T) - dependence follows 
the exponential, i.e., p(T) = po + pi cxp(— Tq/T) rather 
than power-law form; ii) p drops by about an order of 
magnitude when T changes in the range comparable to 
the Fermi energy ej? of 2D electrons; hi) the metallic 
state is quenched by the magnetic field. 

Here, we are not going to discuss existing attempts 
pd| |l8[ to interpret these experiments. (We found in 
Refs. [[ll 18 1 no satisfactory physical explanation of the 
substantial drop in the resistivity in a narrow temper- 
ature interval in an obviously nonsuperconducting sys- 
tem). Instead, we propose a simple mechanism which 
seems to naturally explain all the peculiarities mentioned 
above. We believe that our general idea can be applied to 
all gated semiconductors. However, here we concentrate 
on Si MOS structures, where the important character- 
istics of a 2DEG and of defects are much better known 
than in other systems. 

A typical n-Si MOS structure consists of a metallic 
gate, Si02 layer, and p-type Si substrate. Strong enough, 
positive gate potential attracts electrons which form an 
inversion layer at the Si02/Si interface. It is known Jl~S| ] 
that due to the oxygen deficit in the oxide, there is a 
substantial concentration of defects close to the interface, 
which are capable of trapping charges. Even in state-of- 
the-art devices, there are more than 10 12 hole traps per 
cm 2 , such as Si-Si weak bonds pj| . To introduce the idea 
of our mechanism, we assume all of the hole traps to be ★) 
characterized by the same energy of the electron level St, 



and **) located at the same distance z from the interface. 
We shall abandon assumption **) later on. Effects of a 
finite width of the trap band will be discussed elsewhere. 

At T=0, the trap charge (and spin) state is determined 
by the chemical potential p of the 2DEG. For e t > p, the 
electron level is empty, i.e., a hole is trapped. The trap 
has a charge +e and thus causes strong scattering of 2D 
electrons. It is crucial for our theory that the charge 
state of a trap can be changed by varying the gate volt- 
age V g . Indeed, the bigger V g the smaller e t = e t (V g ). 
At V g = V*{z) determined from e t (V*) = p, the trap 
captures an electron (i.e. emits a hole) and is neu- 
tralized. Being neutral and remote from 2D electrons, 
the defect cannot scatter them any more. Neutraliza- 
tion of the oxide charges reduces resistivity p and thus 
causes an insulator-to-metal transition. When T is high 
(^ \e t — p\), roughly half of the traps are charged. As 
a result, p is rather high and depends weakly on both T 
and V g . On the contrary, for \p, — St\ < T the density 
of charged traps behaves as exp[(e t — p)/T], resulting in 
the exponential /3(r)-dependence [feature (i)]. The tran- 
sition takes place for both degenerate and nondegenerate 
2DEGs [feature (ii)]. Finally, the magnetic field effect 
(iii) can be attributed to the spin freeze-out of holes [^o) : 
Zeeman splitting favors spin 1/2 (charged) state with re- 
spect to the singlet (neutral) state of the defect. 

It should be noted that we neglect here quantum in- 
terference of 2D electrons and thus do not attempt to 
describe the insulating phase. However, we will see that 
even in the classical case dp/dT can change sign due to 
the //(T)-dependence. 

Let us now abandon assumption ★★), i.e., take into ac- 
count a broad distribution of distances z. In order to 
understand why such a distribution does not smear the 
transition, we consider the electrostatic energy of an elec- 
tron in the oxide e e (z). Given the total oxide thickness 
d and its dielectric constant e ox , e e can be written as 



e e {z) = eV g z/d + e 2 /(2e ox z). 



(1) 



Here the two terms represent the external electric field 
and the image force from the 2DEG, respectively (charges 
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induced in the gate can be neglected provided that z <C 



d). e e (z) reaches its maximum e m at z 



where 



e m = -2y/eV g £d\z m = d^J e d /eV g - e d = e /(2e ox d). (2) 

also be expressed through the mean distance r 
between 2D electrons: 



a B r. 



(3) 



where as = r/r s is the effective Bohr radius. Eq.([3j) fol- 
lows from the relation between the 2DEG concentration 
n s = l/(7rf 2 ) and the gate voltage: en s = e ox V g / '(And) 
plf . In order to have a meaning in a macroscopic theory, 
z m has to exceed the screening radius of 2DEG (equal to 
as/4 for Si (001) surface f22|). Therefore in low-density 
devices (r s 3> l/y/2), this length scale is quite legitimate. 

Assuming that the double (hole) occupancy of a trap 
is impossible, the probability of a trap to be charged is 



r 1 

2C 



exp 



' H - e e (z) - e t 
T 



1 



(4) 



where (7 = 1. According to Eqs.(0,^), a homogeneous 
distribution of traps leads to a distribution of charges 
which is peaked at z = z m , the width of the peak being 



6z = d T 2 e d {eV g y 



1/4 



T 2 (eV g e d ) 



1/4 



(5) 



For d — 2000 A and eV g = 1 eV, we get e d — 1 meV and 
e m = 63 meV, so that Sz/z m ~ y/T{K)/l% <C 1, since 
T < 5K. At T = 5 K, z m w 63 A and Sz w 8 A. This 
sharpness of the distribution peak in Eq.([|) manifests 
itself in a sharp metal-insulator transition, as V g is varied. 

How does a positive charge, separated by a distance 
z > from the 2DEG, affect the resistivity? It 
turns out that a bound localized state is formed with 



= *3/4 1/4 



< z and £{, = — e 2 /(e*z) being the local- 



ization length and the energy of this state, respectively 
(e* is the effective dielectric constant of the 2DEG). The 
trap and bound electron form a dipole, which is oriented 
perpendicular to the 2DEG plane. For z ~ as^/rj and 
yfrl ~> 1, the (transport) scattering cross-section E(e, z) 
of such a dipole for electrons with energy e can be eval- 
uated classically: 



E(e,z) = 2.74(e 2 z 2 /2e*e) 1 / 3 . 



(6) 



p= (A t /e 2 n s )V2m*e£(e,z m ) / dzP+(z)( — ) 



The Drude formula for the resistivity can be written as 

rd / Z \ 2 / 3 

dzP+(z)( — 

/o 

(7) 

where iV* is the total volume concentration of the traps, 
m* is the effective mass of the electrons, and e has a 



meaning of their effective energy, which can be expressed 
through the 2DEG Fermi energy e p via 



s F 



^f^) 5/6 cosh- 2 
4T\e F ) 



/ e — ep \ 
V 2T J 



(8) 



Eq. (g) interpolates between two limits: e w 
T/r 6 (ll/6) « 1.44T for T > e F , while in the opposite 
limit e = ef- For T ~ ef, the effective energy e « 2T. 
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FIG. 1. Scaling function 1Z [Eq. @] vs dimensionless tem- 
perature t for several dimensionless gate voltages v [Eq. (M)]. 
u increases in the direction of the arrow. ef/Sci = 0.25. (a) 
Case A. v = -0.2 . . . 0.7. T a /e d = 0.04. Dot-and-dashed line 
indicates the transition, (b) Case B. v — —1.4 . . . 0.7. 

In the saddle-point approximation, Eq.(Q) reduces to 

p=(h/e 2 )pon(V gi T); (9a) 

po = QA6^r- s (e*/2e ox ) 1/6 (N t f/Trn s )(f/d) 2 / 3 ; (9b) 

dx 



/T 3 e\V6 

1{V g ,T) /q [/(r)/2]exp(x 2 + s) + 1 



(9c) 

In Eq. (]9cJ), we took into account the /i(T)-dependence: 
,„ - e t + M(0)1 IT- f(T) = e [MT)-M0)]/T_ (10) 
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We have to consider two distinct cases: (A) chemical 
potentials of the 2DEG and of the Si substrate coincide; 
(B) the 2DEG is disconnected from the substrate. A 
straightforward calculation gives 

Ia{T) = (T/T a ) 3 / 4 ; / B (T) = 1 - exp(- £F /T) (11) 

where T a is determined by the acceptor concentration 
p0| . Although case B is more likely to occur in a real 
device |23| , we shall concentrate mostly on case A which 
exhibits a clear metal-insulator transition even in a clas- 
sical model; see below. 
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FIG. 2. Case A. Scaling function TZ [Eq. @] vs di- 
mensionless gate voltage v for dimensionless temperatures 
t = 0.1... 0.6 [Eq. @]. t increases in the direction of the 
arrow. Inset: a blow-up of the transition region. 



The exponential part of the 7?.(T)-dependence disap- 
pears when s, Eq. (|l0|), vanishes. This happens at 
V g = Vg = [e t - m(0)P /4e£rf and thus V g c defines the 
transition point. The asymptotic behavior of TZ away 
from this point is 



T 3 e 



1/6 



t/ttQ, for < 1; 
ln 1/2 fl ,for n> 1. 



(12) 



where 51 = [2/ f(T)]e~ s . The distance from the transition 
can be measured by S = (K, — V£)/Vg. Provided (5 2 <C 
4T 2 /(e c ;eV r 9 ), variable s in Eq.(|l0|) acquires a scaling form 



s » ^e d eV g c {6/T) = v/t; t = T/e d ; v = 8^eV g c /e d . (13) 

The 7?.(T)-dependence in the scaling region is shown in 
Fig. [I]. For u 3> t, the system is in the "metallic" 
phase characterized by 1Za,b exponentially decreasing 
with t. Due to the /i(T)-dependence, dlZ^/dt changes 
sign at some v slightly bigger than zero, exhibiting thus 
a metal-to-insulator transition. For larger negative v [not 
shown in Fig. 0(a)], the "^(^-dependence saturates. At 



critical behavior for \S\ < 0.01. This is consistent with 
experiments |l}^]. 

In case B, there are two distinct regions: exponentially 
decreasing and ^-independent TZ b , the crossover between 
the two occuring for |«| ~ t. Quantum interference effects 
should result in localization, converting T-independent p 
into an exponentially diverging one; the metal-insulator 
transition in this case will be discussed elsewhere. 

At the transition, TZa ~ 0.1 and TZb ~ l(cf. Fig. 0). 
At the same time, n s ~ 10 11 cm -2 and f ~ 100 A for 
d = 2 • 10 -5 cm and eV g = 1 eV. Estimating N t f ~ 
10 12 cm -2 and r s ~ 10, we obtain for the resistivities at 
the transition p c A ~ O.lh/e 2 and p c B ~ /i/e 2 . These values 
are within the experimentally observed range ||b. 

As Fig. § shows, the transition between the insulating 
and metallic phases in case A is very well-defined, despite 
the fact that TZ does not solely depend on the scaling 
variable v/t. Closer inspection of the transition region 
(Fig. ||, inset) reveals however that the transition occurs 
over a finite range of v rather than at a single point. 
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FIG. 3. Case A. Data collapse in TZ plotted vs t/\v\. Inset: 
"Duality" between "resistivity" p* = TZa /TVa an d u* = 1/p* 
plotted as a function of S [Eq. (13J)]. Solid downward: p*(S). 
Solid upward: cr*(<5). Dashed: cr*(— S). 

Figure || depicts the (approximate) data collapse for 
TZ plotted as a function of t/\v\ = T/T with T = 
\6\y/edeVg. The inset of Fig. || demonstrates the "du- 
ality" feature, i.e., the symmetry between the resistiv- 
ity p in the insulating phase and the conductivity a in 
the metallic one. Experimentally, a similar collapse was 
achieved for To cx |<5| a [in the quantum phase transition 
theory (QPTT), a = vz\. In all of the experiments, ex- 
cept Ref. ||, a is close to 1, i.e. to our prediction. 

An additional insight comes from nonohmic measure- 
ments. In Ref. 0]c, the dependence of p on the source- 
drain voltage Vsd = £L (where £ is the electric field, 
and L is the source-drain distance) was also found to be 
a scaling one p — p(£/£ ) with £ a cx |<5| b . We believe that 
this ^-dependence can be attributed to simple heating. 



3 



Indeed, the effective temperature of electrons T* is deter- 
mined by the energy balance. For strong enough electric 
field, i.e, when T* > T, and for 2D electrons, 



e£JD(T*)T eph (T*) = (tt/V6)T* 



(14) 



Here D(T) is the diffusion constant of electrons at tem- 
perature T, and T ep h(T) oc T~ p is the relaxation time 
of the electron temperature, which we assume to be de- 
termined by electron-phonon scattering (24). One can 
check that if p{T) and D(T) obey a scaling law l/D(T) oc 
p(T) — F(T/Tq), the £ -dependence of the resistivity is 
also a scaling-like: p = G(£/£q), where £ o = oc \8\ aa 
with a = 1 +p/2, and where function G is obtained by 
solving Eq. (14) for given F. If p = 3 (as it is the case for 
good metals), a = 2.5. Experimental value of a = b/a 
HJ is ~ 2.25. This discrepancy can easily be explained 
by p being smaller than 3. On the other hand, QPTT 
predicts a = 1 + z , i.e. a— 2 at z = 1. 

One can check that the strong heating regime is real- 
ized under the conditions of Rcf. [0c, if 



T ep h > 0.1ms/ [£(mV/cm)] 



(15) 



Strong heating of a 2DEG has recently been observed in 
a Si MOSFET device |25| similar in its parameters to 
that used in Refs. fPH H 
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FIG. 4. Case A. Scaling function 1Z vs temperature t for 
Zeeman splittings E z = (0, 0.04, 0.08, 0.12, 0.16, 0.2) x e d (E z 
increases from the bottom to the top curves), v — 0.1. Inset: 
The same as in Fig. [I|a but for Ez = 0.15ed- 

We now turn to the effect of a magnetic field. Consider 
a hole trap, e.g., a Si-Si weak bond |h|. Such a trap can 
find itself in one of the three states with energies Ei\ 

1 = 1,2,3. For i = 1 two electrons occupy the bond. 
This is supposed to be a neutral (Q = 0), singlet (S = 0) 
state. State 2 (3) has one spin down (up) electron on the 
bond. Accordingly, Q = +1 and S = 1/2 for both states 

2 and 3. A magnetic field splits the doublet: E-i — E^) = 
e t ± Ez-, where Ez is the Zeeman splitting. As a result, 



at given T and V g probability P+ to find a trap in a 
Q = +1 state increases with Ez, i.e., with the magnetic 
field: one should substitute C = cosh(Ez/T) instead 
of 1 into Eq.(4). This results in a magnetoresistivity 
demonstrated in Fig. || 
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